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We estimate the amount of temperature-dependent squeezing and entanglement in the coUective 
excitations of trapped Bose-Einstein condensates. We also demonstrate an alternative method of 
temperature measurement for temperatures much less than the critical temperature (T <C Tc)- 
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I. INTRODUCTION 

Since the first experimental realization of Bose- 
Einstein condensation in trapped atomic gases, a number 
of papers have highlighted the parallels between processes 
in atom optics and nonlinear (photon) optics. Well- 
known examples include four-wave mixing [l| and soliton 
dynamics 01, while the squeezing of matter wave fields 
has also been examined H, 0, U- These latter studies 
consider the way in which the nonlinear interactions be- 
tween the atoms can generate entangled atomic beams, 
for example by spin-exchanging collisions of spinor Bose- 
Einstein condensates. 

Entangled states, to which squeezed states are closely 
related, have been extensively studied in the last few 
years because of their potential application in quantum 
information processing and quantum computing. In this 
paper, we discuss squeezing and entanglement in terms 
of the quasiparticle modes of a trapped gas. Quasipar- 
ticle modes represent collective excitations of the Bose- 
Einstein condensate (BEG), and are in fact one of the 
most fundamental features of its dynamics. They are eas- 
ily generated by applying time-dependent perturbations 
to the trapping potential, and were observed experimen- 
tally quite early on; they have since been studied in some 
detail both theoretically and experimentally. 

The low temperature regime in which the majority 
of the atoms are in the condensate has been studied in 
various experiments @, Q • These experiments revealed 
oscillations of the condensate with almost no damping, 
and the results were in good agreement with predictions 
based on the zero temperature Gross-Pitaevskii equation 
(GPE) m . Excitations at higher temperatures were stud- 
ied in later experiments ^ , where large energy shifts and 
rapid damping rates were observed. While the standard 
theory of elementary excitations due to Bogoliubov 
explains most of the observed effects, it cannot account 
for higher-order processes, such as the Beliaev damp- 
ing. And while finite temperature studies using the GPE 
have explained various properties of BECs observed ex- 
perimentally, they cannot fully describe the evolution of 
collective excitations at higher temperatures illj] since 



the spontaneous part of the Beliaev damping of the ex- 
citations is not included. Theoretical techniques have 
recently been developed that take the quasiparticle in- 
teractions into account and experiments to test the 
predictions are planned for the near future. The main 
motivation of the present paper is to analyse these ex- 
periments on the quantal aspects of quasiparticle mode 
evolution in a BEG. 



We note that Beliaev damping may be understood 
in terms of a nonlinear frequency-mixing mechanism in 
which a Bogoliubov quasiparticle of frequency UJ2 inter- 
acts with the ground state (condensate) and generates 
two quasiparticles of frequency uji that divide the ini- 
tial energy equally (Fig. The mechanism is actu- 
ally based on a four-wave process that masquerades as 
a three-wave interaction analogous to optical parametric 
down-conversion (OPDC) in which a photon in a nonlin- 
ear medium splits into two photons of lower energy. 

We calculate the amount of squeezing for a realistic 
system containing 10000 ^^Rb atoms in a trap with spher- 
ical geometry. We also propose an alternative method for 
determining the temperature of a BEG well below the 
critical temperature by observing the variation of the en- 
velope of the collective excitations. This should be more 
accurate than current techniques based on fitting a ther- 
mal Gaussian profile to the atomic cloud. 



The paper is organised as follows: In Sectional we de- 
scribe the higher-order Hamiltonian that accommodates 
processes such as Landau and Beliaev damping in the 
quasiparticle excitations. We derive equations of mo- 
tion for a quasiparticle mode. In Section IIIII we make 
quantitative estimates for the degree of squeezing and 
entanglement in quasiparticle excitations. In addition 
we calculate the expected "damping rate" and coupling 
between modes at various temperatures which can then 
serve as a temperature measurement calibration for ultra- 
cold atoms. A possible way of detecting squeezed states 
in atom optics is discussed briefly in Section [l VI 
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II. HAMILTONIAN IN THE QUASIPARTICLE 
BASIS 

The many-body Hamiltonian for a system of bosons 
with pairwise interactions can be written in the usual 
second quantised formahsm as, 

^ = E ^^'^^^ + ^ E {^j\V\km)alala,a,„, (1) 

ij ijkrn 

where the matrix elements Hij^^ are given by 

Ht^^ j d^vrd^W'^A^)- (2) 

Here, H'^p ~ ~2m^^ ^trap is the single-particle 
Hamiltonian with a confining potential Vtrap, and the 
basis state wave functions are 'ipi{v). {ij\V\km) denotes 
the matrix element for the interaction potential V{v) be- 
tween atoms. The operators a\ and are the creation 
and annihilation operators for mode i that obey the usual 
Bose commutation relations 

[ai, a]] = 5ij, [flj, flj] = [a|, a]] = 0. 

The Hamiltonian |^ is written in a single-particle ba- 
sis where the operator di annihilates a particle from the 

^ For simplicity from now on we drop the explicit notation in the 
operators 



state with wave function ■0i(r). The wave function ■00 (r) 
describes the condensate, while the remaining functions 
form a complete set orthogonal to the condensate. 

Wc follow the standard procedure and rewrite the 
Hamiltonian |^ in powers of dj by replacing the operator 
So with y/N^ = {N—Nex Y^^ where N is the total number 
of atoms in the trap and N^x = X^i^^o ^1^^ number 
of excited atoms. The quadratic part of the Hamiltonian 
can be diagonalised exactly using the standard Bogoli- 
ubov transformation |l2j| . which allows us to write the 
Hamiltonian in a quasiparticle basis. The quasiparticle 
operators Pi are defined by^ 

where the matrices Uij and Vij , which depend on temper- 
ature and the trap geometry, are calculated numerically 
using the Bogoliubov-de Gennes (BdG) formalism 13]. 

The non-quadratic terms, which affect the energy and 
shape of the condensate, are expected to be small and 
hence dealt with perturbatively. We can, therefore, cal- 
culate the quasiparticle shifts and widths at zero and 
finite temperature. It has been shown that the ef- 
fective Hamiltonian can be written in terms of the quasi- 
particle operator /3i as 



H' = const 



E(^ 

i#0 



h.c. 



(4) 



in which only processes that generally conserve energy shift from first order perturbation theory. The coeffi- 
have been included. The constant term in the equation cients Qjk are determined by 
simply defines the zero of energy and Ae^ is the energy 



mnq^O 

+ {mn\V\qO) [u*„UjqVkrn + U*^VjnUkq + vlqVjnVkm] ■ 



(5) 



where the indices i,j and k are labels that denote the 
quasiparticle energy levels. For a 3-D condensate, the 
index i stands for the quantum numbers n, I and m. 
We note that the coefficients are temperature dependent 
through Uij and Vij . 

The Hamiltonian given by equation Q contains terms 



beyond the standard Bogoliubov approximation, thus we 
are using a fuller description of the BEG that involves in- 
teracting quasiparticles and can take important processes 
such as Landau and Bcliaev damping into account. Be- 
liaev processes occur at zero temperature and are dom- 
inant in the low-temperature regime. Landau processes 
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on the other hand predominate at higher temperatures; 
these processes, in which two quasiparticles colhde to 
form a single quasiparticle cannot occur at zero temper- 



where 

= (^p + Aep) /h, 

It is well known that homogeneous systems have a con- 
tinuous spectrum of excitations; however, for a trapped 
system, the spectrum of states to which the excitation 
can couple is discrete. Furthermore, in the case of a 
spherical trap there is a degeneracy in the azimuthal 
quantum number m, which in turn means that fewer 
resonances are present [13. [l5| . This implies that the 
time evolution of the excitation is dominated by a strong 
coupling to only a few modes. On the other hand, the 
geometry of the trap can be changed by adjusting the fre- 
quencies in the radial and axial directions independently, 
and the selection of a dominant mode in this way has 
been demonstrated In particular, a Beliaev process 
has recently been observed for a scissors mode, where 
one mode is resonantly coupled to two modes of half the 
original frequency [l7|; this is the kind of process that 
we seek to model in the present paper. The equations of 
motion for modes p = 1, 2 in a Beliaev process reduce to 

Pi = -iuOiPi - iV2lP2P\, (8) 

Pi = \uiPI+\v2iPIPi, (9) 
P2 = -iw2/32-i^/3i/3i, (10) 

Pi = iu:2pl+i'^P\pl (11) 

where we have used the fact that an — 1^21/2, and we 
have chosen the phases so that the coefficients are real. 

Equations of the form JHJ-ldJ have been studied in 
quantum optics and we apply them here to study the 
squeezing in the quasiparticle excitations and its tem- 
perature dependence. We note that damping of quasi- 
particle excitations has earlier been described in terms 
of nonlinear mixing of quasiparticle modes p^ : however, 
the present work is quite distinct as the operator nature 
of the quasiparticle annihilation operator Pi is retained, 
and the calculation is not restricted to the quadratic ap- 
proximation. We are therefore extending the work of 



ature because there are no excited quasiparticles. 

From the Hamiltonian Q), the Heisenberg equations 
of motion for Pp and are 

(6) 
(7) 



Ref. |l^, so that spontaneous quantum processes are 
included. 

To give an idea of the order of magnitude of the num- 
bers involved, the coefficient 1^21 equals 1.9 x 10~^ in 
trap units for the following parameters: T = 20 nK, 
N = 10000, iotrap = 27r X 100 Hz, for a spherical trap ge- 
ometry and using ^^Rb atoms for which the s-wave scat- 
tering length a — llOoo, where ao is the Bohr radius. The 
values of 1/21 are significantly dependent on temperature; 
which implies that the equations of motion for Pi and p\ 
are altered accordingly. We can exploit this feature as 
a method of measuring temperature below Tc- Figure [3 
indicates the dependence of 1/21 on temperature, which 
is an effective dependence on the number of particles in 
the condensate. We notice that at higher temperatures 
the coefficient decreases swiftly, whilst for lower temper- 
atures it is fairly constant. 

III. TEMPERATURE DEPENDENT COUPLING 
PROCESS 

A. Non-depleted regime 

We consider first the case where the higher (p = 2) 
mode - the pump - has a much larger population than 
the lower one and where the pump depletion is ignored 
so that the operator P2 can be approximated by a c- 
number 62. Physically, this would represent a situation 
in which the higher mode is being continuously driven by 
a resonant excitation. The solutions of the equations of 
motion for Pi and p\ in the interaction picture are then 
given by 

( \ _ ( cosh(m) -isinh(f7i) \ ( /3i(0) \ 
V Pli^) /V isinh(17i) cosh{m) ) \ p\{{)) ) ' 

(12) 

where — 1^21^2, and we have chosen the phases such 
that the coefficients are real. This approximation ne- 
glects the depletion of the "pump" and the solution will 
cease to be valid once appreciable down-conversion oc- 
curs. 

Solutions (|12|l are formally equivalent to equations that 
describe degenerate parametric down conversion in quan- 



iPp = UlpPp+ ^ (TjkP]Pk+ ^ VjkPlPj, 

j.k^a j.k^a 
-iPl = cjppl + ^ a*^p]pl + ^ ,^*kp]pk, 
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turn optics: a process known to be a significant source of 
squeezed states. In our case tlie squeezing parameter is 
given by 



r = nt, 

which is temperature dependent through O 
Defining the quadrature operators c 



(13) 



and Y{t) 
{AY{t)f 



m = 

we have the variances 



I) 



4 

{2Ni + 1; 



exp(-2m), (14) 
exp(2m), (15) 



where Ni is the number of particles in mode 1. In optics, 
the quadrature operators are well-defined quantities cor- 
responding to the amplitude and phase of the electromag- 
netic (EM) oscillation. One may think of the amplitude 
and (temporal) phase of oscillations in a quasiparticle ex- 
citation, in a similar way, and we therefore interpret our 
quadrature operators X and Y as amplitude and phase 
of oscillations. Equations H14I) and H15|l imply that the 
degree of squeezing depends on the amount of the lower- 
mode atoms present. This is directly related to the tem- 
perature of the ultracold atoms in the condensate [l^ . 
defined in terms of the initial Bose- Einstein (BE) distri- 
bution of quasiparticles. Surface plots of AX and AF are 
shown in Figures 13 a) and b). At t = the behaviour of 
the variances is determined by the BE distribution func- 
tion for mode 1. At higher temperatures, the behaviour 
of the coefficient V21 dominates the trend. 

Experimentally, what is of interest is correlation func- 
tions such as (PiPi) or (/3|/3i), among others. These 
quantities have a number of significant physical interpre- 
tations. Assuming an initial number state in the quasi- 
particle basis, some important correlation functions are 
given by: 

{f3i{t)/3i{t)) = -i(Ni + ^\smh{2nt), (16) 



cosh(2m) - -. (17) 



A surface plot of {(3i(3i) as a function of temperature, 
T and time, t is given in Figure0]a); a non-zero value of 



this quantity implies the presence of squeezing and hence 
entanglement. The correlation function of equation l|17|) 
is just the population in the lower mode, and is directly 
measurable as the amplitude of oscillation. A surface 
plot of this quantity is shown in Figure^b) as a function 
of temperature (T) and time {t). We note that at t — 
the trend of both correlation functions is given by the 
BE distribution function for the population in the lower 
mode. As time evolves, the population of the lower mode 
increases. The solution is valid within the non-depleted 
regime, where mode 2 is continually replenished. 

B. Depleted regime 

We now examine the case where mode 2 is depleted and 
an explicit analytical solution is possible. The analysis 
is handled more easily by considering the equations of 
motion for the number operators Ni = pjPi, which from 
equations (j51- lll|l are 



"dT 

dN2 
dt 



(18) 
(19) 



Equations (|18|l and l|19|) can, in fact, be uncoupled by 
calculating the second derivative and using the fact that 
the operator A — Ni+ 2N2 is a constant of motion. The 
uncoupled equations are given by 



d^iVi 

dt2 



dt2 



i^iii-SN^ + 2AN1+ A), 
-|i(12iV|-8A7V2-FA2- 



A). 



(20) 
(21) 



In the case in which N2 is strong, we can determine the 
population by taking the average of equation (|^ . We 
end up with a c-number second-order differential equa- 
tion with the following initial conditions 



dN2 



dt 



t=o 
N2{0) 

whose solution is given by 



= 0, 
= N20; 



(22) 
(23) 



N20 + ia2 - iV2o)sn2 f-^^-i^t, ^^^) , for N20 < a2, 
N2{t)^{ (24) 

ai + (iV20 - ai)nd2 (m^^I^t, ^ , for N2, > as- 



Here sn(u, fc) and nd(u, fc) are Jacobi elliptic functions |20l l2l|. ai and a2 are the roots of the quadratic poly 
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nomial PiN^) = -iV| + [A - N2o)N2 - A^|o " " 
A) /A + N20A. As we have pointed out, the couphng co- 
efficient 1^21 depends on temperature; consequently the 
behaviour of A'2 must change with T. The initial value 
A^20 was calculated for each temperature according to the 
BE distribution function and an initial driving was also 
considered. These parameters were then used in the solu- 
tion (|24|l . FigureElis a plot of the population A'2 showing 
its dependence on temperature and time. 

In a practical case, the fitting of an experiment to 
the theory can be realized by noticing that the initial 
part of the curve for the population in the quasiparticle 
mode 2 could be approximated by a function of the form 
ai cos(rt) (Fig. EJ. To a first approximation the quantity 
r is given by 

/(ai - A^2o)(A2o - a2) 
'^ = "^^V • ^''^ 

The full solution, given by equation H24II . has been used 
in a fitting routine and a plot of F as a function of tem- 
perature is shown in Figure In principle, F has the 
advantage of being readily measurable. At lower tem- 
peratures the effect of the initial population is dominant, 
while, for higher temperatures, the coupling coefficient 
1^21 decreases rapidly so that the effect of the initial popu- 
lation is masked, which happens at approximately 0.6Tc. 
This effect is the result of having a competition between 
the changing coupling strength and the influence of the 
initial population. It is noted that one could also consider 
F as an experimentally accessible indicator of the amount 
of squeezing in the quasiparticles, as it is proportional to 
the coupling strength 1/21- 

IV. DISCUSSION 

Squeezed states of the electromagnetic field were real- 
ized some time ago and the degree of squeezing can be 



measured with standard techniques such as homodyne 
detection. In the case of squeezed quasiparticle modes 
of a BEC, it will be necessary to devise a method for 
their detection. It is noted that a related work that uses 
neutron scattering was suggested by Yurke . 

We suggest for our system that Raman scattering be- 
tween the two excited phonon states could be used to 
characterise the amount of squeezing present since sin- 
gle particle transitions will be strongly modified by the 
presence of a correlated pair excitation. This effect has 
been demonstrated for a homogeneous gas by Stamper- 
Kurn et.al. |23j | who used a single particle transition to 
excite a phonon. We then have two routes for the trans- 
fer of momentum. The first moves a particle with zero 
momentum to momentum fc, while the second takes the 
momentum from —k to zero. The interference between 
the two routes will infiuence the single particle transi- 
tion. Raman scattering will be strongly infiuenced by 
the presence of coherent states in the lowest mode. The 
precise form will depend on the conditions, pulses used, 
and other experimental parameters p^ . Similar ideas 
have been used to probe coherent phonons in solid state 
physics psj . 

In summary, we have analysed a process analogous 
to parametric down-conversion in trapped Bose-Einstein 
condensates. Physically this corresponds to Beliaev 
damping, where two quanta are produced by one of 
higher frequency. This, in turn, suggests the possibil- 
ity of observing temperature-dependent squeezing in the 
collective excitations of the condensate. The amount of 
squeezing gives a direct indication of the temperature and 
vice versa. 
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Figure Captions 

1. Schematic diagram of energy levels for a trapped 
Bose-Einstein condensate. In the Beliaev process, 
analogous to optical parametric down-conversion, 
a quasiparticle of frequency uj2 interacts with the 
ground state generating two quasiparticles of fre- 
quency oJi that divide the initial energy equally. 



2. Coefficient 1^21 vs temperature for TV = 10000, 
^trap = 2TT X 100 Hz, a spherical trap geometry 
and using ^^Rb atoms. 

3. a) Surface plot of AX as a function of temperature 
and time, b) Surface plot of AY as a function of 
temperature and time. At t = 0, the behaviour is 
fully determined by the Bose-Einstein (BE) distri- 
bution function. 

4. Correlation functions calculated within the non- 
depleted regime: a) The average (/3i/3i). The oc- 
currence of squeezing is implied by a non-zero value 
of this quantity, b) Surface plot of (/3|/3i), which 
describes the behaviour of the population in mode 
1. 

5. Surface plot of the population of mode 2 as a func- 
tion of time and temperature. 

6. Plot of N2{t), the solid line is the solution for the 
following parameters: A^20 = 151.4225, at T = 
20nK. The dashed line represents the amplitude 
of oscillations of the quasiparticle mode 2. 

7. The parameter F plotted against temperature is ob- 
tained so that the curve ai cos(rt) corresponds to 
the envelope shown in figure IHI 
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